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Abstract. We find statistically significant spatial correlations between the arrival 
directions of the highest energy cosmic rays (HECRs) observed by the Akeno Giant Air 
Shower Array (AGASA) and large-scale structure of galaxies observed by Sloan Digital 
Sky Survey (SDSS) in the redshift ranges of 0.006 < z < 0.012 and 0.012 < z < 0.018 
at angular scale within ~ 5°. This result supports a hypothesis that the sources of 
HECRs are related to galaxy distribution even in the northern sky, which has been 
already indicated by Pierre Auger Observatory in the southern sky. We also investigate 
the dependency of the correlation on the absolute magnitude, color, and morphology of 
the galaxies. For galaxies with 0.006 < z < 0.012, the correlation tends to be stronger 
for luminous and red galaxies. Based on these results, we discuss plausible HECR 
sources and constraint on Galactic magnetic field. 
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1. Introduction 

Although the sources of the highest energy cosmic rays (HECRs) have been poorly 
known, recent observations begin to unveil a part of the nature of HECR origin. The 
first indication to HECR sources was reported by the Akeno Giant Air Shower Array 
(AGASA). The AG ASA reported the anisotropic distribution of cosmic rays (CRs) 
above 4 x 10 19 eV at small angular scale comparable with their angular uncertainty 
to determine arrival directions of primary CRs [lj. The anisotropic feature has been 
interpreted as the existence of point-like extragalactic sources. An important step to 
their origin is the positional correlation between the arrival directions of HECRs with 
energies above ~ 6 x 10 19 eV and nearby astrophysical objects reported by the Pierre 
Auger Observatory (PAO) [21 [3]. These results strongly suggested that HECRs should 
be of extragalactic origin. Several groups have also analyzed the PAO data by using 
catalogs of different astrophysical objects and then have confirmed the correlation with 
matter (or galaxy) distribution of local Universe [H G2, El [7] . Although the same analysis 
of new PAO data as Refs. [2j [3] shows the decreases of the significance level of the 
correlation [8], the correlation with the large-scale structure is still positive [9]. 

The PAO observes HECRs from the southern sky, since it is located in Argentina. 
If HECR sources are extragalactic objects as suggested by the PAO, the correlation 
of HECRs and galaxy distribution is also expected in the northern sky. However, 
such correlation has never been clearly reported, though the correlation of HECRs 
with the supergalactic plane was reported in the northern sky [TOl [TT] . We tested the 
correlation of the AGASA data with galaxies of Infrared Astronomical Satellite Point 
Source Redshift Survey (IRAS PSCz) [12J, but a significantly positive correlation could 
not be found [TJ. There are two reasons why we have not found such correlation when 
we assume that the correlation exists in the northern sky. One is the small exposure of 
the AGASA relative to the PAO. In the period of Refs. [21 E], the total exposure of the 
PAO is ~ 7 times as large as that of the AGASA. Thus, the AGASA may not have the 
number of highest energy (HE) events large enough to find correlation. The other is the 
completeness of a galaxy catalog used in [7J. In other words, it is the possibility that 
the IRAS PSCz catalog of galaxies does not reflect the large-scale structure of galaxy 
distribution in local Universe. 

In this study, we focus the latter reason, because the former problem cannot be 
solved by us without accumulating more number of HECR events. We investigate the 
correlation between HECR events observed in the northern sky and galaxy distribution 
which has higher completeness than the IRAS PSCz catalog. For this purpose, we 
adopt galaxies observed by Sloan Digital Sky Survey (SDSS) [13]. The SDSS is a survey 
project of galaxies to investigate the large-scale structure of the Universe in optical 
bands and therefore can observe faint galaxies down to m r ~ 20.0 where m r is the 
r-band apparent magnitude of galaxies. Thus, SDSS galaxies better trace the structure 
of the local Universe. As HECR events observed in the northern sky, we adopt the 
AGASA data published in Ref . [H] . 
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This paper is laid out as follows: in Section [21 we explain the data samples of 
HECR events and galaxies in detail. Also, a statistical method which we adopt to test 
the correlation between them is described. In Section [3], we calculate cross-correlation 
functions defined in Section [2] between the AGASA events and SDSS galaxies with 
several redshift ranges. We also investigate the dependencies of the correlation on the 
r-band absolute magnitude, color, and morphology of galaxies. These give us useful 
information to search for the HECR sources because the features of the host galaxies 
reflect the nature of objects contained in the hosts. Section @] is devoted to discuss and 
interpret the results obtained in Section [31 Finally, we conclude this study in Section [51 

2. Data Samples and Statistical Method 

2.1. Sample of Highest Energy Cosmic Rays 

We adopt the AGASA data published in Ref. [14J as a sample of HECR events observed 
in the northern hemisphere. We have two reasons why we use this sample: (i) that the 
AGASA data is the largest sample in the northern sky in which the energy and arrival 
direction of each event are published and (ii) the AGASA is a ground array whose 
directional dependence of the aperture is easier to estimate analytically than that of 
fluorescence telescopes like High Resolution Fly's Eye (HiRes). 

That event sample consists of 57 events with energies above 4 x 10 19 eV in the 
energy-scale of the AGASA, 8 events out of which are so-called super- GZK events 
[T5] which mean that their energies are larger than the energy of a spectral cutoff 
(~ 10 20 eV) predicted from interactions with cosmic microwave background (CMB) 
photons (Greisen-Zatsepin-Kuz'min (GZK) cutoff; [TH Ej). Since these super-GZK 
events are difficult to be simply produced by astrophysical scenarios, their origin have 
often been proposed in beyond standard model physics (for a review, see Ref. [IE]). An 
astrophysical solution is the existence of their sources within several tens of Mpc which 
is comparable with the energy-loss length of photopion production of HECRs in CMB. 
Although the GZK steepening of HECR spectrum is not observed by the HiRes [19] and 
PAO [2D] and is a problem itself, we simply assume that all events are of astrophysical 
origin in this study to use this sample and analyze them. 

2.2. Sample of Galaxies 

We use data from the 7th data release (DR7) of the SDSS [21], especially large-scale 
structure subsample named dr72bright0 sample of the New York University Value Added 
Catalog (NYU-VAGC) [22j. This is a spectroscopic sample of galaxies with u,g,r,i,z- 
band (K-corrected) absolute magnitudes, r-band apparent magnitude m r , redshifts, and 
information on the mask of the survey. Spectroscopic galaxies are selected by conditions 
of 10.0 < m r < 17.6 and 0.001 < z < 0.5. Strictly speaking, the upper bound of 
apparent magnitude of galaxies is less than 17.6 in a very small fraction of regions, but 
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Figure 1. Rcdshift vs. r-band absolute magnitude of SDSS galaxies within z = 0.03. 



we have practically no problem to regard the upper bound as 17.6 because the fraction 
of such regions is less than 0.1%. 

We construct an approximately volume-limited sample of SDSS galaxies to analyze 
the cross-correlation between the AGASA events and galaxy distribution. Fig. [1] plots 
the r-band absolute magnitude of all the SDSS galaxies in the sample within z < 0.03 as 
a function of their redshifts. We can clearly see a line due to the flux limit determined 
by m r = 17.6. 

We study the cross-correlation with galaxies within z = 0.024 which corresponds to 
~ 100 Mpc in the concordance cosmology. The GZK horizon of protons with the energy 
of ~ 4 x 10 19 eV are larger than ~ 100 Mpc. Fig. [2] shows the cumulative fraction of the 
flux of protons arriving from sources with distances representing the horizontal axis to 
the total flux of protons. The left and right panels are for protons with energies of 10 19,6 
eV and 10 19 ' 8 eV, respectively. The cumulative fractions are calculated under several 
assumptions on the nature of HECR sources: a spectral index a and a cosmological 
evolution factor of HECR source density with a shape of (1 + z) m , the static emission of 
protons, and neglecting intergalactic magnetic field (IGMF). For these calculations, we 
assume a ACDM cosmology with H = 71 km s -1 Mpc -1 , Q m = 0.3, and = 0.7, as 
the Hubble parameter, the ratio of matter density to the critical density, and the ratio of 
cosmological constant to the critical density, respectively. We can see that the fraction 
of protons arriving at the Earth from sources within 100 Mpc is only 10-20 %. Even 
for protons with 10 19 ' 8 eV, the fraction is maximally ~ 40 %. However, IGMF could 
more deflect the trajectories of HECRs emitted from sources with larger distances and 
therefore the cross-correlation signal between arriving HECRs and such sources may be 
weaken. In the case of the PAO data, we cannot find any significant signal of positive 
correlation with a galaxy distribution with z > 0.018, whereas significantly positive 
correlation was found for z < 0.018 [7]. The angular scale of the positive correlation is 
up to ~ 15°. These indicate a possibility that HECRs emitted from distant sources are 
significantly deflected by intervening magnetic fields. Therefore, it is natural to focus on 
the cross-correlation of UHECRs with nearby galaxy distribution in order to constrain 
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Figure 2. Cumulative ratios for sources to contribute to the fluxs of extragalactic 
protons with 10 196 eV (left) and 10 19 8 eV (right) at the Earth. 



their sources. 

In order to construct a volume-limited sample of galaxies, we consider galaxies with 
the absolute magnitude below -17 within z = 0.024. Under this criterion, all the galaxies 
with M r < —17 within 100 Mpc are included in the galaxy sample, as we can see in Fig. 

m 

2.3. Cross-correlation Function 

We adopt a cross-correlation function used in our previous work [7] , which was originally 
introduced in Ref. [23] for the cross-correlation between two types of galaxies, 

EG(9)-EG'(9)-E'G(9) + E'G'(9) 

= . (1) 

E, G, E' and G' represent HECR events, galaxies in the catalog, HECR events randomly 
put with number density proportional to the detector aperture, and galaxies randomly 
put following the angular selection function (e.g., survey window, bright star mask, and 
so on) in the observed sky, respectively. EG (6) is the number of pairs between E and 
G with the angular distance from 9 to 9 + A9 normalized by the number of HECR 
events multiplied by the number of galaxies. We set A9 = 1°, which is comparable 
to the typical angular determination scale of the arrival directions of primary cosmic 
rays for HECR observatories. EG' (9), E'G(9), and E'G'(9) are also defined similar to 
EG(9). E' and G' enable us to correct the non- uniformities of the HECR aperture and 
the galaxy sampling. 

The aperture of a ground array depends on the declination of observed directions 
reflecting the daily rotation of the Earth. The declination dependence of the exposure 
(= aperture x observation time) can be analytically estimated as [21] 

ui(5) cx cos(a ) cos(<5) sin(a m ) + a m sin(a ) sin(5), (2) 
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where a m is given by 
and 



if £ > 1 

7T if £ < -1 

cos _1 (£) otherwise 



(3) 



cos(6>) — sin(a ) sm(S) 



(4) 



cos(ao) cos(<5) 



Here, oq is the terrestrial latitude of a ground array and 6 is the zenith angle for a data 
cut because of experimental reasons. For the AGASA, a = 35°47' and 9 = 45° [TJ]. We 
produced the E' sample in Eq. (1) following Eq. (2). We set the number of E' sample 
to be twenty times more than the actual HECR events. 

The survey area of the SDSS is also limited and very complicated due to masks. The 
SDSS observes about one fifth of the whole sky and the depth of observations slightly 
depends on directions. We take account of these effects by distributing random galaxies 
(i.e., G' in Eq. 1) following the angular selection function. We use mangle software 
[25| [26] to do so. Again we distribute at least twenty times more random galaxies than 
the real ones, which is an usual choice, and then confirmed that this particular choice 
does not affect the results. 

3. Results 

This section is dedicated to show results of our analysis. First of all, we consider the 
cross-correlation between the AGASA events and SDSS galaxies in 4 ranges of redshifts 
within z = 0.024, which corresponds to ~ 100 Mpc in the concordance cosmology, in 
Section 13.11 This analysis will lead to a result that the AGASA events correlate more 
with SDSS galaxies in intermediate redshift ranges. Next, we focus on these galaxies 
and investigate the dependency of the cross-correlation on the absolute magnitude of 
the galaxies 13.21 Then, we also consider the dependencies of color and morphology of 
galaxies to find the nature of environment for HECR generation 13.31 

3.1. Dependence on Redshift 

We divide our sample of galaxies into 4 parts in redshift for every 0.006, which 
corresponds to ~ 25 Mpc. Note that the flux of HECRs from sources in each shell 
is comparable if we assume that HECR composition is proton and IGMF is neglectable 
as we can see in Fig. [2], i.e., all the lines in Fig. [2] can be well approximated by straight 
lines within 100 Mpc. We consider the cross-correlations between the AGASA data and 
these 4 samples of galaxies. Results are shown in Fig. [31 

In the left panels, we show the positions of galaxies {red dots) and the arrival 
directions of the AGASA events with E > 4 x 10 19 eV (circles) for 4 redshift ranges 
of galaxies: 0.001 < z < 0.006, 0.006 < z < 0.012, 0.012 < z < 0.018, and 
0.018 < z < 0.024 from top to bottom, for visibility. The super-GZK events, which are 
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Figure 3. left: Skymaps of the AGASA events and SDSS galaxies with M r < —17. 
The redshift ranges of SDSS galaxies (plotted in red) considered are < z < 0.006, 
0.006 < z < 0.012, 0.012 < z < 0.018, and 0.018 < z < 0.024 from top to bottom. 
The AGASA events are also ploted as circles (E > 4 x 10 19 eV; green, E > 10 20 
eV; blue), right: Cross-correlation functions between the AGASA events and SDSS 
galaxies calculated from the corresponding left panels {histogram). We also show the 
cross-correlation functions calculated from randomly distributed events in green. 
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Table 1. Numbers of galaxies used in calculations of cross-correlation functions. 



events with energies beyond 10 20 eV, are shown as blue circles. The numbers of galaxies 
we used are summarized in Table [H We can see that the distribution of galaxies is 
highly structured. For instance, Virgo cluster is seen at the center of the top panel. 
Any AGASA event does not spatially correlate with the positions of Virgo cluster, as 
seen in our previous works jTj, which is the same situation as the PAO events [27] . Also, 
we can find Coma cluster at the center of the bottom left panel and the AGASA event 
neither correlates with Coma cluster. However, we see that many event correlates with 
filamentary structures of galaxies by eye. This may be quantified by the cross-correlation 
function. 

The right panels of Fig. E are the cross-correlation functions between the AGASA 
events and the positions of galaxies shown in the corresponding left panels (histograms) . 
In order to compare them to the cases of isotropic events, we also show the cross- 
correlation functions between isotropically distributed events and the galaxies in green. 
We simulate isotropically distributed events with the same number as the AGASA data 
100 times and estimate the averages (crosses) and standard deviations of the cross- 
correlation functions (error bars). These standard deviations can be regarded as 1 a 
statistical errors. 

In the top panel, there is no signal indicating a positive correlation and the signal 
is consistent with isotropic events, which reflects the fact that the AGASA events do 
not correlate with not only Virgo cluster but also with the other structures in the 
nearest Universe. Since the deflections of HECRs from nearby sources are expected 
to be small, there might not be HECR sources in the Universe within ~ 25 Mpc 
in the northern sky. In intermediate redshift ranges, which are 0.006 < z < 0.012, 
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Figure 4. Dependence of the cross-correlation on the r-band absolute magnitude of 
SDSS galaxies with the redshift of 0.006-0.012. 



0.012 < z < 0.018, we can find positive signals within ~ 5°. These reflect about 
10 events seen to correlate with filamentary structures of galaxies. These correlations 
are similar to the PAO results in Ref. [2j [3] because a large fraction of the PAO- 
correlated AGNs, which are a representative of large-scale structure of local Universe, 
is distributed at around z ~ 0.01, which corresponds to ~ 50 Mpc [28]. Also, local 
Universe with such a size may have a volume enough to contain HECR sources and the 
deflections of UHECR trajectories may not be critical. Thus, we can interpret these 
positive signals as the appearance of HECR sources. Under this interpretation, the 
angular scale of the positive correlations corresponds to the maximum value of HECR 
deflections in the local Universe. Finally, positive correlation is not found in the bottom 
panel. Large deflections of HECRs due to their large propagation distance might let 
correlation disappear. 

3.2. Dependence on Absolute Magnitude 

In the previous section, we found the positive correlation between the AGASA events 
and SDSS galaxies within 0.006 < z < 0.012 and 0.012 < z < 0.018. In this section, 
we focus on galaxies in these redshift ranges and investigate the dependence of the 
cross-correlation signals on their r-band absolute magnitude. 
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Figure 5. Same as Fig. [4] but the redshift range of galaxies is 0.012 < z < 0.018. 



Fig. H] shows the cross-correlation functions of the AG ASA data and SDSS galaxies 
with 0.006 < z < 0.012 in several ranges of their r-band absolute magnitude. Although 
SDSS galaxies with M r < —17.0 are taken into account for a volume-limited sample 
in Section [37TI we also consider galaxies with —17.0 < M r < —15.0 here (upper panel) 
to investigate the correlation with faint galaxies. Note that M r = —15.0 is smaller 
than the flux limited by m r = 17.6 and therefore these galaxies are completed in this 
redshift range (see Fig. [1]). The lower panels show the cross-correlation functions with 
galaxies with —19.0 < M r < —17.0 and M r < —19.0 as histograms. We also show the 
cross-correlation functions with the galaxies and isotropic HECR distribution with 57 
events, and their standard deviations (~ la statistical errors), as for Fig. [3] We can 
find the strongest positive signal at the angular scale within ~ 5° in the 3 panels in the 
lower right panel. Although the number of galaxies in this M r range is only 145, we can 
see such strong signals against a simulation by isotropic events. This implies that the 
AGASA events were generated in the luminous (generally large) galaxies or at least in 
an environment where there is a number of luminous galaxies. 

Fig. [5] shows the same as Fig. HI but the redshift range considered is 0.012 < 
z < 0.018. Since the r-band absolute magnitude corresponding to the limited flux is 
smaller than in 0.006 < z < 0.012, we do not consider an absolute magnitude range of 
— 17.0 < M r < —15.0. Instead of this, we consider galaxies with —17.0 < M r < —16.0 
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Figure 6. Numbers of SDSS galaxies with the redshifts of 0.006 < z < 0.012 (left) 
and 0.012 < z < 0.018 (right) as a function of a concentration parameter, Ci nv , and the 
difference of absolute magnitude at (/-band and r-band. The numbers increase with 
approaching red from blue. 



for fainter galaxies. We cannot find clear difference like Fig. [4] among the 3 ranges of 
absolute magnitude. Reasons might be a smaller flux of HECRs for each source and 
larger deflection angles of their propagating trajectories due to source distances larger 
than the case of Fig. HI In order to investigate whether the same trend as in Fig. H] is 
realized in this redshift range or not, we need more number of HECR events. 

3.3. Dependence on Color and Morphology 

The color of a galaxy reflects the population of stars in the galaxy and therefore the 
evolution history of the galaxy. Thus, the dependence of HECR-galaxy correlation on 
the color of galaxies informs us of the epochs of HECR generation in the history of 
galaxies and environment around HECR sources. This can be a hint to constrain the 
sources of HECRs. 

Colors of galaxies are generally defined as difference between magnitudes in two 
different bands in astronomy, and therefore the definition of colors is not unique. Here, 
we define the colors of galaxies as M g — M r , which is difference between (/-band and 
r-band absolute magnitudes. Since the wavelength of the r-band is larger than g-band, 
a galaxy with large M g — M r is called relatively red and vice versa. M g — M r is a good 
indicator of contents of galaxies because the 4000A break of galaxies is included in g- 
band. This break is the appearance of long-lived stars like K-type stars. Thus, strong 
star forming and a number of young (massive) stars are not expected in galaxies with a 
strong spectral break at ~ 4000A, i.e. relatively red galaxies for the color definition. 

The morphologies of galaxies also have information on objects in the galaxies. 
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Figure 7. Dependence of the cross-correlation on the color of SDSS galaxies with the 
redshift of 0.006-0.012. 



Galaxies accompanying strong non-thermal radio emission, like radio galaxies, have 
a tendency to be early- type (elliptical) morphology [29|. On the other hand, several 
works have pointed out that the host galaxies of 7-ray bursts (GRBs) might have late- 
type (like irregular or spiral) morphology j30j[3T]. The NYU-VAGC attaches the ratio 
of two Petrosian radii c- mv = r^o/rgo called a (inverse) concentration index. This index is 
a good indicator to the morphological type of each galaxy [32]. Since elliptical galaxies 
tend to fade out slowly with the distance from the center, they have relatively smaller 
concentration indices. On the other hand, late-type galaxies have large concentration 
indices. Ref. [32] compared concentration indice with morphologies of SDSS galaxies 
with their apparent magnitude above 16.0 detected during a commissioning phase 
and found that the best criterion between early-type and late-type morphologies is 
Cinv ~ 0.35. 

Fig. [6] shows the frequency distributions of galaxies with M r < — 17 and the redshift 
ranges of 0.006 < z < 0.012 (left) and 0.012 < z < 0.018 (right) on a parameter space 
of Ci n v and M g — M r . The number of galaxies increases when the color in each cell is 
close to red. In both panels, the distribution of galaxies can be separated above and 
below M g — M r ~ 0.6. Thus, we define galaxies with M g — M r above and below 0.6 as 
red and blue galaxies in this study, respectively. On the other hand, the distribution 
cannot be clearly separated in the direction parallel to c- mv , but a large fraction of the 
galaxies have c inv > 0.35, which means late-type galaxies. 

Fig. [7]shows the cross-correlation functions of the AGASA data with SDSS galaxies 
with 0.006 < z < 0.012 which have M g - M r > 0.6 (red; left) and M g - M r < 0.6 
(blue; right). We can obviously find that relatively red galaxies strongly correlate with 
the AGASA events, while the correlation with blue galaxies is consistent with random 
simulations. This implies that HECRs are generated in red galaxies or at least in an 
environment in which many relatively red galaxies exist. Since a red galaxy is luminous 
and thus large because it is after many stars are made, this result is consistent with that 
in the previous section. Such a feature is not clear for galaxies with 0.012 < z < 0.018 
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Figure 8. Same as Fig. but the redshift range of galaxies is 0.012 < z < 0.018. 



as seen in Fig. [HI 

We specially focus on galaxies with M r < — 19 and 0.006 < z < 0.012, which 
strongly correlate with the AGASA events. Fig. [9] shows plots of the galaxies on 
a parameter space of c imr and M g — M r (red). We also plot galaxies within 5° from 
any AGASA events in green, whose angular scale reflects the result that the spatial 
correlation scale is ~ 5° (see Fig. [3]). We can find that all the galaxies within 5° from 
the AGASA events are M g — M r > 0.6, relatively red, while c im are broadly distributed. 
We also check the morphologies of these galaxies (green points) by eye by using SDSS 
DR7 navigate tool [33] . 2 galaxies with c imr < 0.35 are elliptical, but the other have 
late-type morphologies. 

In the summary of this section, we found that the AGASA events significantly 
correlate with SDSS galaxies with 0.006 < z < 0.012 and 0.012 < z < 0.018 within the 
angular scale of ~ 5°. Especially in the case of 0.006 < z < 0.012, the correlation 
tended to be stronger for luminous (M r < —19) and red galaxies. The luminous 
galaxies around the AGASA events were relatively red and morphologically late- type. 
For 0.012 < z < 0.018, the dependence of the correlation on colors and morphologies 
of galaxies was not clear. This was possible because the flux of HECRs from a sources 
with this redshift range is smaller than for 0.006 < z < 0.012 and/or the trajectories of 
HECRs are more deflected than for sources with 0.006 < z < 0.012. These dependencies 
should be tested in the future when the more number of HECR events are accumulated. 



4. Discussion 



In this section, we discuss the sources of HECRs based on results found in Section O We 
specially focus on the results on galaxies with 0.006 < z < 0.012, though it is non-trivial 
that the correlation features in 0.006 < z < 0.012 are universal. Potential sources of 
HECRs ever proposed are GRBs [3U [35], [36], newly born magnetars (37J, radio-loud 
active galactic nuclei (AGNs) [38] [39] HQ] , and clusters of galaxies [U][l2]. We discuss 
the relations between our results and the features of the host galaxies of GRBs (Section 
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Figure 9. Distribution of galaxies with 0.006 < z < 0.012 and M r < —19 on a Ci nv - 
M g — M r plane. The red crosses are all the galaxies and the green circles are galaxies 
which are placed within 5° from the AG AS A events. 



14. ip . magnetars (Section l4.2p . and AGNs (Section I4.3P under the assumption that the 
observed correlation is true correlation between HECRs and their sources. The case of 
clusters of galaxies is not discussed. In Section 14.41 we refer to the possibility that the 
correlation is fake due to the propagation of HECRs in magnetized intergalactic space 
for completeness. Finally, we discuss a constraint on Galactic Magnetic Field (GMF) 
from the results in this study in Section 14.51 

4-1. j-ray Bursts 

Since GRBs (at least long-duration ones) are thought to be related to explosions of 
massive stars at the end of their life [I3j 0U US], their host galaxies are expected to have 
high star formation rate. Although the GRB population with redshift is small, many 
observations have indicated that the typical nature of GRB host galaxies is, in fact, of 
faint star-forming galaxies dominated by a young stellar population [46] . In general, the 
host galaxies have low luminosities [30], low masses |31j . and therefore low metallicities 

HUSH]. 

Although the host galaxies are too faint to be observed by galaxy surveys, it is 
no problem in local Universe where we consider because the absolute magnitudes of 
galaxies observed in the local Universe are comparable with those of GRB host galaxies 
observed in high redshift Universe. Fig. 6 of Ref. [31 J gave the histogram of the number 
of GRB host galaxies as a function of B-band (centered at 442 nm) absolute magnitude. 
Almost all of GRB host galaxies have —22 < Mb < —17. This has been unchanged in 
a recent larger database [H] Since a large fraction of SDSS galaxies has M g — M r > 0, 
GRB host galaxies are expected to be more luminous in r-band than B-band. Note that 
the central wavelengths of g-band and r-band are 477.0 nm close to B-band and 623.1 
nm, respectively. Thus, our approximately volume-limited sample of galaxies, which 
consists of galaxies with M r < —17, includes GRB host candidates. We can directly 
compare our results with the features of GRB host galaxies. 
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The color of GRB host galaxies tends to be blue [21] even in nearby Universe 
[50] . We should notice the definition of the color of galaxies. In Ref. [31], the authors 
discussed the colors of GRB host galaxies by using infrared wavelengths, R — K, which is 
different from the definition in this study. However, the essence of both definitions is the 
same, i.e. these are indicators of old (long-lived) star population, and therefore we can 
compare results of this study with those of GRB host studies qualitatively. It might be 
difficult that GRBs are HECR sources from the qualitative discussion, since the AGASA 
events correlate with relatively red galaxies. Of course, qualitative discussion should be 
performed under the same definition of galaxy colors. 

Since GRBs are transient phenomena, we must consider the time-delay of HECRs 
in intergalactic space and it should be compared with the typical lifetime of massive 
stars as progenitors of GRBs. If the typical time-delay of HECRs is larger than the 
lifetime of massive stars, galaxies which correlate with the AGASA events do not seem 
to be blue. Although the masses of GRB progenitors have uncertainty, the masses are 
larger than ~ 25M Q or typically ~ 40M Q under a collapser model [51] in which requires 
black hole formation by stellar core collapse. A typical main sequence lifetime is ~ 10 7 
yrs. Several IGMF models predict that a fraction of HECRs has time-delay larger than 
~ 10 7 yrs even for protons (52], [53] However, the trajectories of such HECRs are highly 
deflected and the correlation between HECRs and the positions of their sources is lost. 
The correlation with matter distribution can be observed, but it is fake correlation with 
their true sources [M] (see Section WM in detail). In this case, it is difficult to confirm 
GRB origin of HECRs from information on colors of the correlated galaxies. Note that 
other IGMF models predict much smaller time-delay [551 EE]. If such a small IGMF 
is realized in the Universe, the discussion on colors of the correlated galaxies above is 
allowed. 

GRB host galaxies are morphologically irregular or late- type [31]. In our study, Fig. 
[6] shows that a large fraction of galaxies has c\ m > 0.35, which means they are late-type 
galaxies. The completeness of the division is ~ 90% [32]. Thus, GRB origin of HECRs 
is not inconsistent with the results of the analyses. Note that a recent study which use 
the largest sample of GRB host galaxies reports no compelling evidence that GRB host 
galaxies are peculiar galaxies in the point of their morphologies [49]. Thus, unfortunatly 
it is not clear from the results of this study whether GRBs are HECR sources or not. 

4-2. Magnetars 

Magnetars are thought to be neutron stars with extremely strong magnetic fields 
(~ 10 15-16 G), but their production mechanism is still controversial (see [57] for a 
review). Although we know about ten of magnetar candidates, the positions of a large 
fraction of magnetars are associated with supernova remnants or massive stellar cluster 
[58] . This indicates that magnetars are produced through core-collapse of massive stars 
like GRBs. Observationally, GRB060218 (sometimes categolized into X-ray flush) was 
a possible candidate of a GRB leaving a magnetar [59] and its progentor mass was 
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estimated as ~ 2OM , relatively smaller than the masses of typical GRB progenitors. 
Also, it was theoretically proposed that a class of GRBs, called low-luminosity GRBs, 
accompany the production of magnetars [60J. 

The nature of the host galaxies of magnetars is poorly known because almost all 
the magenetar candidates are located in our Galaxy. One of the 2 exceptions is in 
Large Magellanic Cloud (SGR 0526-66) [61] and the other is in Small Magellanic Cloud 
(CXOU J0100-72) [62]. Therefore, there is little information on host galaxies other than 
our Galaxy. However, based on a hypothesis that magnetars are created by core-collapse 
of massive stars as discussed above, they are more frequently produced in star-forming 
galaxies, similar to the case of GRBs, but a little redder galaxies are allowed because the 
masses of magnetar progenitors are thought to be smaller than those of GRBs. Thus, 
the constraints may match the possible conditions of the host galaxies of magnetars a 
bit better than for GRBs. 

4-3. Active galactic nuclei 

Many models to accelerate particles to the highest energies at AGNs have been proposed: 
at giant radio lobes of Fanaroff- Reify II (FRII) galaxies with Mpc-scale jets [38J, 
giant flares of blazars [40], and the shocks of cocoon expansion of small version (~ 1 
kpc) of radio-loud galaxies [39]. An important point is that all of these scenarios 
requires galaxies dominated by non-thermal radiation, which have generally early-type 
morphologies [29], though their required sizes of sources are different. On the other 
hand, a large fraction of SDSS luminous galaxies with 0.006 < z < 0.012 has late- 
type morphologies. Thus, it might be difficult that these source candidates are HECR 
sources. 

4-4- Propagation effects 

In the 3 previous sections, we treat the cross-correlation between HECRs and their 
sources as true one. However, the arrival directions of HECRs can be far from the 
positions of HECR sources without reducing the correlation with matter distribution, 
if IGMF is strong. Note that IGMF strength and structures are still under debate. 

Ref. |63j pointed out the possibility of fake correlation. In this scenario, even the 
HECRs are scattered by strongly magnetized structure of the Universe called scattering 
centers like clusters of galaxies, large magnetized lobes of radio galaxies and so on. 
Therefore, if there is a real source of HECRs far from a region which observed HECRs 
correlate with, the HECRs seems to come from the region when the HECRs are scattered 
by a strong magnetic field in the region. Ref. [63] reported that about 50% of the PAO 
events published in Ref. [3] could be such fake correlation. A similar phenomena is 
observed in Ref. [54] by simulations of propagation of protons considering a sophisticated 
and simulated IGMF model [53]. They also show that predicted angles between the 
arrival directions of HE protons and AGN-like objects defined in their simulations are 
consistent with the correlation angles derived by the PAO [3] even if there are real 
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sources far from the AGN-like objects. In any case, regions which correlate with the HE 
events could not always point out real sources. We should keep it in mind that there are 
much uncertainty on IGMF modelling. There are also IGMF models to predict smaller 
deflections of HECRs [551 [56]. 

4-5. Constraint on GMF 

The angular scale of the positive correlation enables us to constrain intervening 
magnetic fields, especially Galactic magnetic field (GMF). Since the structure of GMF 
is independent of galaxy distribution in local Universe, GMF can reduce the expected 
spatial correlation between the arrival directions of HECRs and the positions of nearby 
galaxies if GMF largely deflects HECR trajectories. However, we found the correlation 
of HECRs and large-scale structure of matter distribution. Thus, GMF is not too large 
to weaken such correlation. 

Moreover, we can constrain GMF structure from the results of this study under 
the assumption that the observed correlation is true correlation between HECRs 
and their sources, because the spatial correlation is sensitive to GMF models in 
the northern sky [6^[ [65]. Refs. [6H [65] investigated the correlation between HE 
protons and their sources taking their deflections by several GMF models into account. 
The angular scale of positive correlation which we found in Section [3] was estimated 
as ~ 5°. This is comparable with the angular scale of positive correlation with 
HECR sources expected by bisymmetric spiral GMF models predicted in those works. 
Although axisymmetric GMF models predicted cross- correlation functions consistent 
with isotropically distributed HE protons within error bars, this does not reject the 
possibility of the axisymmetric GMF models. However, protons with energies of 10 19,8 
eV are deflected by more than 5° in about a half of the northern sky for the axisymmetric 
GMF models. Thus, it might be difficult that the axisymmetric models are realized, 
i.e., there might be reversals of GMF at least outside the solar system. Furthermore, 
the small deflection implies proton-dominated composition of HECRs. 

We should notice 2 limitations on the discussion above. One is that the SDSS 
surveys about a half of the northern sky. The discussion above is based on the cross- 
correlation between HECRs and galaxies in a sky where the SDSS observed. The other 
is the definition of the angular scale for positive correlation in this study. That angular 
scale is determined as angular scale in which correlation signals exceed to the 1 o 
error bars attached into the cross-correlation function calculated from random event 
distribution. When the number of observed events increases, the error bars decrease 
and the angular scale of the positive correlation might be larger. 

5. Conclusion 

In this study, we investigated the correlation between the arrival directions of HECRs 
detected in the northern sky (the AGASA data) and galaxy distribution observed by 
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the SDSS. The SDSS can observe faint galaxies, and therefore SDSS galaxies are a 
good sample to investigate the correlation with the large-scale structure of matter in 
local Universe. We found the positive correlation of the AGASA events with SDSS 
galaxies within the angular scale of ~ 5° with the redshift ranges of 0.006 < z < 0.012 
and 0.012 < z < 0.018 for the first time. This result supported a hypothesis that 
the sources of HECRs are related to galaxy distribution even in the northern sky, 
which has been already indicated by several analyses of the PAO data in the southern 
sky [21 [31 HI El El [7j. We also tested the dependencies of the correlation signals on 
the r-band absolute magnitude, color, and morphology of the galaxies, since these 
information is generally useful to constrain HECR sources involved in the host galaxies. 
For galaxies with 0.006 < z < 0.012, we found that the AGASA events correlated with 
luminous (M r < —19) and red (M g — M r > 0.6) galaxies, and the luminous galaxies 
around the AGASA events were relatively red and morphologically late-type. Such 
dependencies were not clear for galaxies with 0.012 < z < 0.018. This was possible 
because the flux of HECRs from a sources with this redshift range is smaller than for 
0.006 < z < 0.012 and/or the trajectories of HECRs are more deflected than for sources 
with 0.006 < z < 0.012. 

We discussed possible source candidates based on our results for galaxies with 
0.006 < z < 0.012 in Section HI Unfortunately, the interpretation of the results 
is ambiguous at present because of uncertainty of relations between HECR source 
candidates and the nature of their host galaxies, and our poor knowledge on intervening 
magnetic fields. Under the limitations, possible interpretation is summarized as follows: 
Morphology dependence of galaxies correlating with the AGASA events is positive for 
GRBs and magnetars, but the colors of these galaxies are relatively red, contrary 
to relatively blue colors expected for these objects. Note that the discussion on 
the colors was only qualitative. More detail and quantitative discussion is required 
to confirm/reject GRB or magnetar origin of HECRs. On the other hand, HECR 
generation by non-thermal phenomena related to AGNs might be difficult from results 
in this study because the morphologies of luminous galaxies near the AGASA events 
are relatively late-type, though such phonomena are strongly related to galaxies with 
early-type galaxies. 

The interpretation above is based on the assumption that the trajectories of HECRs 
are not largely deflected by IGMF. When IGMF is strong, fake correlation could occur 
[631 EI]- The correlation between the arrival directions of HECRs and the positions of 
their sources is reduced, though the correlation with matter distribution is conserved. In 
this case, it is not possible to confirm HECR sources by the nature of galaxies correlating 
with HECR events. 

The angular scale of the positive correlation enabled us to constrain intervening 
magnetic fields under the assumption that we see the correlation between HE protons 
and their sources. The angular scale estimated in this study, < 5°, was compared with 
theoretical predictions by Ref. [65]. The angular scale was consistent with predictions 
based on bisymmetric spiral field models of GMF. Since the larger angular deflections of 
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HE protons are predicted for axisymmetric GMF models, it might be difficult to realize 
the axisymmertic models. 

In order to investigate correlation signals between HECRs and matter 
distribution/their sources more accurately, we need more number of HECR events. It 
was not clear whether the correlation features found in this study for galaxies with 
0.006 < z < 0.012 are common in the Universe. This can be checked by increasing the 
observed number of HECRs. We can also understand intervening magnetic fields and the 
nature of HECR sources. Ref. [66J predicted the spatial correlation between HE protons 
and their sources in local Universe even taking a realistically structured IGMF model 
into account, which was clearly visible by accumulating ~ 200 events above ~ 6 x 10 19 
eV for the number density of HECR sources of 10~ 5 Mpc -3 . Data accumulation in the 
northern sky is in progress. The total exposure of the Telescope Array have already 
reached 75% of the AGASA exposure [67J and is expected to reach that of the AGASA 
in this year. Thus, their data is quite useful to test the correlation and the nature 
of the host galaxies of HECR sources. Also, two projected HECR observatories with 
extremely large effective area, Extreme Universe Space Telescope (JEM-EUSO) [68J and 
the northern site of the PAO [69], enable us not only to clarify the correlation but also 
to unveil the positions of HECR sources in the northern sky of local Universe. 
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